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ABSTRACT 

The  seasonal  ocean  circulation  and  the  seasonal  thermal  structure  in  the  South  China  Sea  (SCS)  were  studied 
numerically  using  the  Princeton  Ocean  Model  (POM)  with  20-km  horizontal  resolution  and  23  sigma  levels 
conforming  to  a  realistic  bottom  topography.  A  16-month  control  run  was  performed  using  climatological  monthly 
mean  wind  stresses,  restoring-type  surface  salt  and  heat,  and  observational  oceanic  inflow/outflow  at  the  open 
boundaries.  The  seasonally  averaged  effects  of  isolated  forcing  terms  are  presented  and  analyzed  from  the 
following  experiments:  1)  nonlinear  dynamic  effects  removed,  2)  wind  effects  removed,  and  3)  open  boundary 
inflow/outflow  set  to  zero.  This  procedure  allowed  analysis  of  the  contribution  of  individual  parameters  to  the 
general  hydrology  and  specific  features  of  the  SCS:  for  example,  coastal  jets,  mesoscale  topographic  gyres,  and 
countercurrents.  The  results  show  that  the  POM  model  has  the  capability  of  simulating  seasonal  variations  of 
the  SCS  circulation  and  thermohaline  structure.  The  simulated  SCS  surface  circulation  is  generally  anticyclonic 
(cyclonic)  during  the  summer  (winter)  monsoon  period  with  a  strong  western  boundary  current,  a  mean  maximum 
speed  of  0.5  m  s_1  (0.95  m  s_1)>  a  mean  volume  transport  of  5.5  Sv  (10.6  Sv)  (Sv  =  106  m3  s_1)>  and  extending 
to  a  depth  of  around  200  m  (500  m).  During  summer,  the  western  boundary  current  splits  and  partially  leaves 
the  coast;  the  bifurcation  point  is  at  14°N  in  May  and  shifts  south  to  10°N  in  July.  A  mesoscale  eddy  on  the 
Sunda  shelf  (Natuna  Island  eddy)  was  also  simulated.  This  eddy  is  cyclonic  (anticyclonic)  with  maximum  swirl 
velocity  of  0.6  m  s-1  at  the  peak  of  the  winter  (summer)  monsoon.  The  simulated  thermohaline  structure  for 
summer  and  winter  are  nearly  horizontal  from  east  to  west  except  at  the  coastal  regions.  Coastal  upwelling  and 
downwelling  are  also  simulated:  localized  lifting  (descending)  of  the  isotherms  and  isohalines  during  summer 
(winter)  at  the  west  boundary.  The  simulation  is  reasonable  when  compared  to  the  observations.  Sensitivity 
experiments  were  designed  to  investigate  the  driving  mechanisms.  Nonlinearity  is  shown  to  be  important  to  the 
transport  of  baroclinic  eddy  features,  but  otherwise  insignificant.  Transport  from  lateral  boundaries  is  of  con¬ 
siderable  importance  to  summer  circulation  and  thermal  structure,  with  lesser  effect  on  winter  monsoon  hydrology. 
In  general,  seasonal  circulation  patterns  and  upwelling  phenomena  are  determined  and  forced  by  the  wind,  while 
the  lateral  boundary  forcing  plays  a  secondary  role  in  determining  the  magnitude  of  the  circulation  velocities. 


1.  Introduction 

The  South  China  Sea  (SCS)  is  a  semienclosed  tropical 
sea  located  between  the  Asian  landmass  to  the  north 
and  west,  the  Philippine  Islands  to  the  east,  Borneo  to 
the  southeast,  and  Indonesia  to  the  south  (Fig.  1),  with 
a  total  area  of  3.5  X  106  km2.  It  includes  the  shallow 
Gulf  of  Thailand  and  connections  to  the  East  China  Sea 
(through  the  Taiwan  Strait),  the  Pacific  Ocean  (through 
the  Luzon  Strait),  the  Sulu  Sea  (through  the  Mindoro 
Strait),  the  Java  Sea  (through  the  Gasper  and  Karimata 
Straits),  and  to  the  Indian  Ocean  (through  the  Strait  of 
Malacca).  All  of  these  straits  are  shallow  except  the 
Luzon  Strait  whose  maximum  depth  is  2400  m.  Con- 
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sequently  the  SCS  is  considered  a  semienclosed  water 
body  (Huang  et  al.  1994).  The  complex  topography  in¬ 
cludes  the  broad  shallows  of  the  Sunda  shelf  in  the 
south/southwest;  the  continental  shelf  of  the  Asian  land- 
mass  in  the  north,  extending  from  the  Gulf  of  Tonkin 
to  Taiwan  Strait;  a  deep,  elliptical  shaped  basin  in  the 
center ;  and  numerous  reef  islands  and  underwater  pla¬ 
teaus  scattered  throughout.  The  shelf  that  extends  from 
the  Gulf  of  Tonkin  to  the  Taiwan  Strait  is  consistently 
near  70  m  deep,  and  averages  150  km  in  width;  the 
central  deep  basin  is  1900  km  along  its  major  axis 
(northeast-southwest)  and  approximately  1100  km 
along  its  minor  axis,  and  extends  to  over  4000  m  deep. 
The  Sunda  shelf  is  the  submerged  connection  between 
southeast  Asia,  Malaysia,  Sumatra,  Java,  and  Borneo 
and  is  100  m  deep  in  the  middle;  the  center  of  the  Gulf 
of  Thailand  is  about  70  m  deep. 

The  SCS  is  subjected  to  a  seasonal  monsoon  system 
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Fig.  1.  Geography  and  isobaths  showing  the  bathymetry  (m)  of  the  South  China  Sea. 


(Wyrtki  1961).  From  April  to  August,  the  weaker  south¬ 
westerly  summer  monsoon  winds  result  in  a  wind  stress 
of  over  0.1  N  m~2  (Fig.  2a)  that  drives  a  northward 
coastal  jet  off  Vietnam  and  anticy clonic  circulation  in 
the  SCS  (Fig.  3a).  From  November  to  March,  the  stron¬ 
ger  northeasterly  winter  monsoon  winds  corresponds  to 
a  maximum  wind  stress  of  nearly  0.3  N  nrr2  (Fig.  2b) 
causing  a  southward  coastal  jet  and  cyclonic  circulation 
in  the  SCS  (Fig.  3b).  The  transitional  periods  are  marked 
by  highly  variable  winds  and  surface  currents. 

The  observed  circulation  patterns  of  the  intermediate 
to  upper  layers  of  the  SCS  are  primarily  forced  by  the 
local  monsoon  systems  (Wyrtki  1961),  with  contribu¬ 
tions  from  the  Kuroshio  Current  via  Luzon  Strait  (or 
called  Bashi  Channel),  in  the  southern  half  of  the  Luzon 
Strait.  The  Kuroshio  enters  the  SCS  through  the  south¬ 
ern  side  of  the  channel,  then  executes  a  tight,  anticy- 
clonic  turn  and  exits  the  SCS  near  Taiwan.  The  volume 
transport  through  Luzon  Strait  is  uncertain.  Different 
authors  reported  values  ranging  from  8-10  Sv  (Sv  = 
106  m3  s-1)  (Huang  et  al.  1994)  to  3  Sv  (Wyrtki  1961). 
This  flow  may  exert  a  strong  influence  on  the  properties 


of  the  northern  SCS  waters.  Hu  and  Liu  (1992)  believed 
this  flow  to  contribute  to  currents  in  the  Taiwan  Strait; 
however,  Li  et  al.  (1996)  showed  that  this  flow  did  not 
contribute  to  currents  in  the  Taiwan  Strait  using  a  nu¬ 
merical  model. 

Eddy  behavior  in  the  SCS  has  two  distinct  features. 
First,  the  number  of  cold  eddies  is  far  greater  than  the 
number  of  warm  eddies  and,  second,  the  eddies  are  sig¬ 
nificantly  affected  by  the  bottom  topography  and  are 
most  likely  to  occur  near  localized  high  current  velocity 
(Huang  et  al.  1994).  Small-scale  eddies  with  seasonal 
dependence  have  been  found  off  coastal  Vietnam  in 
summer,  near  Natuna  Island  and  in  reef  areas.  Large- 
scale  eddies  have  been  found  primarily  during  the  sum¬ 
mer  monsoon  (Huang  et  al.  1994).  Based  on  more  com¬ 
plete  datasets  such  as  the  U.S.  Navy’s  Master  Ocean¬ 
ographic  Observational  Dataset  and  National  Centers 
for  Environmental  Prediction  monthly  sea  surface  tem¬ 
perature  (SST)  dataset,  Chu  et  al.  (1997a,  b;  1998)  iden¬ 
tified  multieddy  structure  in  the  SCS. 

Upwelling  and  downwelling  occur  off  the  coast  of 
central  Vietnam  and  eastern  Hainan.  The  summer  mon- 
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(a) 


(b) 


Fig.  2.  Climatological  wind  stress  for  (a)  Jun  and  (b)  Dec  (after  Hellerman  and  Rosenstein 

1983). 


soon  is  believed  to  cause  an  Ekman-type  drift  current 
and  corresponding  offshore  transport,  leading  to  up- 
welling  (Huang  et  al.  1994).  Others  (Chen  et  al.  1982) 
have  pointed  out  that  upwelling  also  occurs  in  these 
areas  during  winter  monsoon  conditions,  casting  doubt 
on  previous  ideas  concerning  the  origin  of  these  fea¬ 
tures. 

In  the  north,  the  waters  are  cold  and  saline.  The  annual 
variability  of  salinity  is  small,  due  to  the  inflow  and 
diffusion  of  high  salinity  water  from  the  Pacific  Ocean 
through  Luzon  Strait.  In  the  south  the  tropical  conditions 
cause  the  waters  to  be  warmer  and  fresher.  During  the 
transitions  the  central  region  is  alternately  subjected  to 
high  and  low  salinity  inflow  as  the  monsoons  reverse, 
resulting  in  a  region  of  higher  horizontal  gradient  and 
annual  variability.  Mixed  layer  depths  vary  from  30  to 
40  m  during  the  summer  monsoon,  and  70  to  90  m 
during  the  winter  monsoon  with  the  variation  due  to 
both  wind  and  current  (Wyrtki  1961). 

Three-dimensional  ocean  models  have  been  used  sev¬ 
eral  times  in  the  past  to  simulate  the  circulation  in  the 
SCS.  Integrating  a  semi-implicit  12-layer  shallow-water 
model  with  50-km  resolution  under  both  mean  winter 
and  summer  conditions  for  15  days,  Pohlmann  (1987) 
simulated  the  reversal  of  the  upper-layer  circulation  be¬ 
tween  the  summer  and  winter  monsoon  seasons.  Inte¬ 
grating  a  23-level  primitive  equation  model  developed 
at  Princeton  University  (Blumberg  and  Mellor  1987) 
with  20-km  horizontal  resolution  under  monthly  mean 
climatological  wind  forcing  (Hellerman  and  Rosenstein 
1983)  for  a  year,  Chu  et  al.  (1994)  simulated  the  seasonal 
variation  of  the  SCS  circulation  and  Chu  and  Chang 
(1997)  explained  the  formation  of  the  SCS  warm-core 
eddy  in  boreal  spring.  Integrating  a  global  1.5 -layer  re¬ 


duced  gravity  model  with  Vi  resolution,  Metzger  and 
Hurlburt  (1996)  successfully  simulated  the  upper-layer 
circulation  and  the  mass  exchange  between  the  SCS,  the 
Sulu  Sea,  and  the  Pacific  Ocean. 

The  objective  of  this  study  is  to  simulate  the  SCS 
thermohaline  structure,  as  well  as  the  circulation,  and 
to  investigate  physical  processes  causing  seasonal  vari¬ 
ability.  We  used  the  Princeton  Ocean  Model  (POM)  to 
examine  the  mechanisms  causing  seasonal  variation  of 
the  SCS  circulation  and  thermal  structure.  The  control 
run  is  designed  to  best  simulate  reality  against  which 
each  experiment  is  compared.  In  the  experiments,  var¬ 
ious  external  and  internal  factors  are  modified  and  the 
resulting  circulation  patterns  and  magnitudes  compared 
to  the  control  run  results.  Specifically  we  estimate  the 
contribution  (in  terms  of  volume  transport,  sea  surface 
elevation  and  circulation  patterns)  of  nonlinear  advec- 
tion,  wind  forcing,  and  lateral  boundary  transport  to  the 
ocean  features  identified  in  the  control  results.  From  this 
we  can  estimate  the  relative  importance  of  these  factors 
to  SCS  oceanography. 

2.  The  numerical  ocean  model 

a.  Model  description 

Coastal  oceans  and  semienclosed  seas  are  marked  by 
extremely  high  spatial  and  temporal  variability  that  chal¬ 
lenge  the  existing  predictive  capabilities  of  numerical 
simulations.  The  POM  is  a  time-dependent,  primitive 
equation  circulation  model  on  a  three-dimensional  grid 
that  includes  realistic  topography  and  a  free  surface 
(Blumberg  and  Mellor  1987).  Tidal  forcing  was  not  in¬ 
cluded  in  this  application  of  the  model  since  high  fre- 


2974 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  29 
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Fig.  3.  Observational  surface  circulation:  (a)  Jun  and  (b)  Dec  (after  Wyrtki  1961). 
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quency  variability  of  the  circulation  is  not  considered. 
River  outflow  is  also  not  included.  However,  the  sea¬ 
sonal  variation  in  sea  surface  height,  temperature,  sa¬ 
linity,  circulation,  and  transport  are  well  represented  by 
the  model.  From  a  series  of  numerical  experiments,  the 
qualitative  and  quantitative  effects  of  nonlinearity,  wind 
forcing,  and  lateral  boundary  transport  on  the  SCS  are 
analyzed,  yielding  considerable  insight  into  the  external 
factors  affecting  the  region  oceanography.  The  hori¬ 
zontal  spacing  of  0.179°  by  0.175°  (approximately  20- 
km  resolution)  and  23  vertical  sigma  coordinate  levels. 
The  model  domain  is  from  3.06°S  to  25.07°N,  98.84°E 
to  121.16°E,  which  encompasses  the  SCS  and  the  Gulf 
of  Thailand,  and  uses  realistic  bathymetry  data  from  the 
Naval  Oceanographic  Office  DBDB5  database  (5'  by  5' 
resolution).  As  pointed  out  by  one  of  the  anonymous 
reviewers,  the  DBDB5  data  are  inaccurate  especially  in 
shallow  regions  comparing  to  the  Defense  Mapping 
Agency  (DMA)  maps.  We  will  modify  the  model  to¬ 
pography  files  to  reflect  what  DMA  maps  indicate  in 
future  studies. 

Consequently,  the  model  contains  125  X  162  X  23 
horizontally  fixed  grid  points.  The  horizontal  diffusiv- 
ities  are  modeled  using  the  Smagorinsky  (1963)  form 
with  the  coefficient  chosen  to  be  0.2  for  this  application. 
The  bottom  stress  rb  is  assumed  to  follow  a  quadratic 
law 

T*  =  P0CD\\bWb,  (1) 

where  p0  (=1025  kg  m-3)  is  the  characteristic  density 
of  the  sea  water,  is  the  horizontal  component  of  the 
bottom  velocity,  and  CD  is  the  drag  coefficient  which 
is  specified  as  0.0025  (Blumberg  and  Mellor  1987)  in 
our  model. 

b.  Atmospheric  forcing 

The  atmospheric  forcing  for  the  SCS  application  of 
the  POM  includes  mechanical  and  thermohaline  forcing. 
The  wind  forcing  is  depicted  by 

du  dv\ 

— ,  —  =  (r0x,  T0y),  (2) 

dz  dzL=0 

where  («,  v)  and  (r0x,  r0y )  are  the  two  components  of 
the  water  velocity  and  wind  stress  vectors,  respectively. 
The  wind  stress  at  each  time  step  is  interpolated  from 
monthly  mean  climate  wind  stress  (Hellerman  and  Ro¬ 
sens  tein  1983),  which  was  taken  as  the  value  at  the 
middle  of  the  month.  The  wind  stress  has  a  typical  mag¬ 
nitude  of  0. 1-0.2  N  m-2  (Fig.  2).  Over  the  two  monsoon 
seasons  the  wind  varies  with  location  and  time,  leading 
to  a  complicated  distribution  of  wind  stress. 

Surface  thermal  forcing  is  depicted  by 

30  (Qh\ 

kh-  =  +  *2C(0obs  -  0)  (3) 

dS 

Ks—  =  (*,(),+  <*2C(Sohs  -  S ),  (4) 

dz 


where  6ohs  and  Sohs  are  the  observed  potential  temper¬ 
ature  and  salinity,  Cp  is  the  specific  heat,  and  QH  and 
Qs  are  surface  net  heat  and  salinity  fluxes,  respectively. 
The  relaxation  coefficient  C  is  the  reciprocal  of  the  re¬ 
storing  time  period  for  a  unit  volume  of  water.  The 
parameters  (a,,  a2)  are  (0,  l)-type  switches:  a,  =  1, 
a2  =  0,  would  specify  only  flux  forcing  is  applied;  a, 
=  0,  a2  =  1,  would  specify  that  only  restoring  type 
forcing  is  applied.  In  this  study,  the  surface  thermal 
forcing  is  determined  solely  by  restoring  forcing;  that 
is,  a,  =  0  and  a2  —  1  in  (3)-(4).  The  relaxation  co¬ 
efficient  C  is  taken  to  be  0.7  m/day,  which  is  equivalent 
to  a  relaxation  time  of  43  days  for  an  upper  layer  30  m 
thick  (Chu  et  al.  1996).  The  net  effect  is  to  prevent  any 
deviation  from  climatology  and  ensure  that  the  SCS  acts 
as  a  heat  source.  The  mixing  coefficients  KM ,  KH,  and 
Ks  were  computed  using  a  level  two  turbulence  closure 
hypothesis  (Mellor  and  Yamada  1982). 

c.  Lateral  boundary  forcing  (Wyrtki  type) 

Closed  lateral  boundaries,  that  is,  the  modeled  ocean 
bordered  by  land,  were  defined  using  a  free- slip  con¬ 
dition  for  velocity  and  a  zero  gradient  condition  for 
temperature  and  salinity.  No  advective  or  diffusive  heat, 
salt,  or  velocity  fluxes  occur  through  these  boundaries. 

At  open  boundaries,  the  numerical  grid  ends  but  the 
fluid  motion  is  unrestricted.  Uncertainty  at  open  bound¬ 
aries  makes  marginal  sea  modeling  difficult.  Three  ap¬ 
proaches,  local-type,  inverse-type,  and  nested  basin/ 
coastal  modeling  are  available  for  determining  the  open 
boundary  condition.  Here,  we  take  the  local-type  ap¬ 
proach,  that  is,  to  use  the  radiative  boundary  condition 
with  specified  volume  transport.  When  the  water  flows 
into  the  model  domain,  temperature  and  salinity  at  the 
open  boundary  are  likewise  prescribed  from  the  cli¬ 
matological  data  (Le vitus  1984).  When  water  flows  out 
of  the  domain,  the  radiation  condition  was  applied: 

d  d 

~(e,  S )  +  U—(6,  S)  =  0,  (5) 

at  on 

where  the  subscript  n  is  the  direction  normal  to  the 
boundary. 

As  we  pointed  out  in  the  introduction,  flows  through 
SCS  straits  are  quite  uncertain.  It  is  not  an  easy  job  to 
choose  one  among  various  estimations.  Despite  being 
old,  Wyrtki’s  (1961)  data  provide  a  balanced  estimation 
of  volume  transports  for  the  Luzon  Strait,  Taiwan  Strait, 
and  Gasper-Karimata  Strait  with  seasonal  variations 
(Table  1).  Since  there  are  no  reliable  estimations  at  the 
Balabac  Channel,  Mindoro  Strait,  and  Strait  of  Malacca, 
we  assumed  zero  transport  there.  Such  a  treatment,  es¬ 
pecially  at  the  Mindoro  Strait,  may  distort  the  solution. 

To  overcome  this  model  weakness,  a  logical  approach 
is  to  develop  a  nested  basin/coastal  model.  Lateral  open 
boundary  conditions  will  be  obtained  from  the  basin 
model  output.  For  the  moment,  we  may  call  our  lateral 
open  boundary  forcing  as  Wyrtki- type  forcing. 
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Table  1 .  Bimonthly  variation  of  volume  transport  (Sv)  at  the  lateral 
open  boundaries.  The  positive/negative  values  mean  outflow/inflow 
and  were  taken  from  Wyrtki  (1961). 


Month 

Feb 

Apr 

Jun 

Aug 

Oct 

Dec 

Gaspar-Karimata 

Straits 

4.4 

0.0 

-4.0 

-3.0 

1.0 

4.3 

Luzon  Strait 

-3.5 

0.0 

3.0 

2.5 

-0.6 

-3.4 

Taiwan  Strait 

-0.9 

0.0 

1.0 

0.5 

-0.4 

-0.9 

d.  Initial  conditions  and  initialization 

The  model  was  integrated  with  all  three  components 
of  velocity  ( u ,  v,  w)  initially  set  to  zero,  and  with  tem¬ 
perature  and  salinity  specified  by  interpolating  clima¬ 
tology  data  (Le vitus  1984)  to  each  model  grid  point. 
The  model  year  consists  of  360  days  (30  days  per 
month),  day  361  corresponds  to  1  January.  It  was  found 
that  90  days  were  sufficient  for  the  model  kinetic  energy 
to  reach  quasi-steady  state  under  the  imposed  conditions 
(Fig.  4).  In  order  to  first  capture  the  winter  monsoon  to 
summer  monsoon  transition,  the  model  was  started  from 
day  300  (30  October),  and  run  to  day  390  (30  January 
next  year)  for  spinup.  After  day  390,  the  model  was  run 
another  three  years  for  each  experiment. 

e.  Mode  splitting 

For  computational  efficiency,  the  mode  splitting  tech¬ 
nique  (Blumberg  and  Mellor  1987)  is  applied  with  a 
barotropic  time  step  of  25  seconds,  based  on  the  Cour- 
ant-Friederichs-Levy  (1928)  computational  stability 
(CFL)  condition  and  the  external  wave  speed;  and  a 
baroclinic  time  step  of  900  seconds,  based  on  the  CFL 
condition  and  the  internal  wave  speed. 

f  Experiment  design 

Our  approach  was  to  carry  out  four  numerical  ex¬ 
periments:  one  control  and  three  sensitivity  runs.  All 
runs  were  completed  for  the  same  three  year  period 
encompassing  both  summer  and  winter  monsoons,  and, 
except  as  specified  below,  utilized  the  same  initial  con¬ 
ditions  (on  day  300). 

Run  1  was  the  control  run  (Wyrtki-type  lateral  open 
boundary  forcing).  The  three  sensitivity  runs  were  run 
2:  linear  dynamics,  run  3:  no  winds,  and  run  4:  zero 
lateral  transport  at  the  open  boundaries.  The  difference 
between  the  control  and  the  sensitivity  runs  at  each  grid 
point  and  time  should  then  isolate  the  nonlinear  dynam¬ 
ics,  wind  forcing  effect,  and  lateral  boundary  effect, 
enabling  independent  examination  of  each. 

The  thermohaline  structure,  sea  surface  elevation,  cir¬ 
culation  patterns,  and  volume  transport  that  constitute 
the  SCS  oceanography  will  be  identified  from  the  con¬ 
trol  run  for  subsequent  comparative  analysis.  An  im¬ 
portant  assumption  is  made  that  the  differences  are  lin- 


Julian  Day  relative  to  Jan  1 


Fig.  4.  Temporal  variation  of  total  kinetic  energy.  The  model  is 
integrated  from  1  Nov  and  reaches  quasi-steady  state  after  90  days 
of  integration. 


ear,  that  is,  higher-order  terms  and  interactions  are  neg¬ 
ligible  and  can  be  ignored.  This  assumption  will  be 
examined  and  qualified. 

3.  The  simulation  (control  run) 

a.  Circulation 

The  most  obvious  feature  of  both  the  summer  and 
winter  SCS  circulation,  measured  and  modeled,  is  the 
western  boundary  current,  the  Vietnam  coastal  jet 
(VCJ).  Hinted  at  in  Wyrtki’s  (1961)  depiction  but  more 
explicit  in  our  model  are  the  cross-basin  currents  located 
at  1 1°-14°N.  The  model  simulates  both  the  summer  and 
winter  SCS  circulation  quite  well  compared  to  the  ob¬ 
servational  study  (Wyrtki  1961). 

During  the  summer  monsoon  period  (mid-May  to  Au¬ 
gust)  winds  blow  from  the  southwest  and  the  SCS  sur¬ 
face  circulation  generally  follows  suit  with  anticyclon- 
icity  in  the  southern  basin  (Fig.  5a).  Inflow  is  through 
the  southern  Gaspar  and  Karimata  Straits  and  outflow 
is  through  the  northern  Taiwan  Strait  and  eastern  Luzon 
Strait  (Table  1).  The  simulated  summer  (June- August) 
mean  general  circulation  pattern  has  the  following  fea¬ 
tures.  Velocities  reach  1  m  s_1  at  the  peak  of  the  summer 
monsoon  within  the  VCJ.  The  western  boundary  current 
splits  into  two  currents  at  12°N:  the  coastal  current  and 
the  offshore  current.  The  offshore  current  further  bi¬ 
furcates  and  partially  leaves  the  coast;  the  bifurcation 
point  is  near  the  south  Vietnam  bight  at  10°N,  110°E. 
The  cross-basin  zonal  current  reaches  14°N  and  a  core 
speed  of  0.3  m  s_1.  The  coastal  branch  continues  north, 
then  east  at  Hainan  Island. 

During  the  winter  monsoon  period  (November  to 
March)  the  winter  Asian  high  pressure  system  brings 
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(a)  (b)  (c) 


Longitude  (E)  Longitude  (E)  Longitude  (E) 


Fig.  5.  Surface  circulation  for  control  run  on  (a)  15  May,  (b)  15  Jul,  (c)  15  Sep,  (d)  15  Oct,  (e)  15  Dec,  and  (c)  15  Mar. 


strong  winds  from  the  northeast  and  the  SCS  surface 
circulation  pattern  is  cyclonic  (Figs.  5d-f).  Inflow  from 
Luzon  Strait  (the  Kuroshio  intrusion)  and  from  the  Strait 
of  Taiwan  augments  currents  southwestward  along  the 
Asian  continental  shelf,  then  southward  along  the  coast 
of  Vietnam  and  eventually  out  through  the  Gasper  and 
Karimata  Straits  in  the  south.  Western  intensification  of 
the  general  cyclonic  circulation  pattern  was  also  sim¬ 
ulated.  From  the  south  coast  of  Hainan  Island,  the  cur¬ 
rent  intensifies  as  it  flows  from  north  to  south  along  the 
Vietnam  coast.  Average  speed  is  around  0.8  m  s_1  in 


the  core.  As  winter  progresses,  near-surface  currents 
along  the  coast  of  Borneo  begin  to  turn  northwest,  even¬ 
tually  flowing  directly  away  from  the  coastline  and  into 
the  Natuna  eddy  (at  5°N,  110°E)  where  the  southern 
edge  of  the  deep  basin  meets  the  Sunda  shelf  (Fig.  5e); 
the  south  extension  of  the  VCJ  veers  westward.  In  early 
spring  a  northward  surface  current  developed  adjacent 
to  the  VCJ  (Fig.  5f)  and  persists  through  the  transition 
periods  to  the  start  of  the  summer  monsoon. 

The  magnitude  of  the  Ekman  surface  current,  VE0 ,  is 
calculated  by 
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Fig.  6.  Mean  surface  circulation  for  control  run  during  (a)  summer  and  (b)  winter. 


\Z1ttt 

EQ  DEp0 1212  sin^l’ 

where  r  is  the  magnitude  of  the  surface  wind  stress,  12 
the  angular  speed  of  rotation  of  the  Earth,  <p  the  latitude, 
and  De  the  Ekman  layer  depth.  For  a  rough  estimation 
on  the  effects  of  Ekman  drift  on  the  SCS  surface  cur¬ 
rents,  we  use  the  following  values: 

0  =  12°,  De  =  100  m. 

The  maximum  wind  stress  t  in  the  winter  (summer)  is 
0.3  N  m-2  (0.1  N  m-2).  The  maximum  surface  Ekman 
current  is  around  0.43  m  s-1  (0.14  m  s-1)  in  the  winter 
(summer). 

In  order  to  compare  the  difference  of  the  two  mon¬ 
soon  seasons,  we  averaged  the  model  variables  such  as 
velocity,  surface  elevation,  temperature,  and  salinity 
over  the  two  periods:  June- August  (summer)  and  De- 
cember-February  (winter.)  The  surface  circulations  for 
the  summer  (Fig.  6a)  show  the  strong  VCJ  flowing 
northward  to  northeastward.  The  winter  circulations 
(Fig.  6b)  indicate  a  reverse  pattern:  the  VCJ  flowing 
southward  to  southeastward.  The  seasonal  variation  of 
VCJ  is  illustrated  by  the  zonal  cross  sections  of  v  com¬ 
ponent  at  13°N  (Fig.  7)  for  winter  (December-February) 
and  summer  (June- August).  The  summer  VCJ  is  flow¬ 
ing  northward  with  a  mean  maximum  speed  of  0.5  m 
s_1  and  a  width  of  100  km,  extending  to  a  depth  of 
around  200  m  (Fig.  7a).  The  winter  VCJ  is  a  much 
stronger  southward- flowing  boundary  current  with  a 
mean  maximum  velocity  of  0.95  m  s-1  and  a  width  of 
100  km,  extending  to  a  depth  of  500  m  (Fig.  7b). 

Average  summer  sea  surface  elevation  varies  from 
—  0.1  to  0.1  m  with  maximum  values  at  the  northeastern 
part  near  Fuzon  (Fig.  8a.).  Sea  surface  heights  for  the 


winter  show  a  0.2-m  depression  over  the  deep  basin 
with  maximum  values  in  the  southwestern  part  over  the 
Sunda  shelf  and  the  Gulf  of  Thailand  (Fig.  8b).  These 
results  are  consistent  with  the  harmonic  analysis  on  TO- 
PEX/Poseidon  altimetry  data  (Yanagi  et  al.  1997).  Fig¬ 
ure  9  shows  the  cotidal  and  corange  charts  of  the  annual 
tide  constituent,  Sa.  The  amplitude  of  Sa  is  large  at  the 
northeastern  and  southwestern  parts  (Fig.  9a)  and  its 
phase  differs  180°  on  both  sides  (Fig.  9b).  The  Sa  con¬ 
stituent  takes  its  maximum  in  August  at  the  northeastern 
part  and  in  December  at  the  southwestern  part.  Thus, 
both  numerical  simulation  and  satellite  data  analysis 
suggest  that  such  a  surface  elevation  anomaly  is  caused 
by  the  current  axis  shift  of  the  seasonal  circulation  and 
the  geostrophic  relation. 

b.  Thermohaline  structure 

Isotherms  and  isohalines  for  summer  and  winter  are 
nearly  horizontal  from  east  to  west  except  at  the  coastal 
regions  (Fig.  10).  Focalized  lifting  of  the  isotherms  and 
isohalines  during  summer  at  the  west  boundary  indicates 
coastal  up  welling  inside  the  VCJ  (Figs.  10a  and  10c). 
Focalized  descent  of  the  isotherms  and  isohalines  during 
winter  at  the  west  boundary  indicates  coastal  down- 
welling  inside  the  VCJ  (Figs.  10b  and  lOd).  Both  coastal 
up  welling  and  down  welling,  combined  with  the  high 
velocity  shear  across  the  VCJ,  results  in  baroclinic  in¬ 
stability.  This  mechanism  may  contribute  to  the  summer 
bifurcation  of  the  VCJ. 

Fatitudinal  thermohaline  variation  for  summer  and 
winter  is  also  shown  in  Figs.  lOe-h.  In  general,  model 
thermohaline  structure  is  consistent  with  the  two  SCS 
water  masses  described  by  Ma  (1998.)  Over  the  southern 
basin  there  is  a  general  lifting  of  isotherms  of  40-50  m 
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Fig.  7.  Mean  latitudinal  velocity  (ms  l)  for  control  run  at  the  13°N  cross  section  during  (a) 
winter  and  (b)  summer.  Contour  interval  is  0.1  m  s-1. 


(a)  (b) 


Longitude  (E)  Longitude  (E) 

Fig.  8.  Mean  surface  elevation  (m)  for  control  run  during  (a)  summer  and  (b)  winter.  Contour 

interval  is  0.025  m. 
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Sa  Constituent  Sa  Constituent 


Fig.  9.  Corange  and  cotidal  charts  of  Sa  constituent  (after  Yanagi  et  al.  1997). 


from  winter  to  summer  above  200  m.  In  the  northern 
SCS,  near-surface  waters  (20-50  m)  are  influenced  by 
the  winter  inflow  of  North  Pacific  Kuroshio  water.  They 
are  uniformly  colder  and  more  saline,  and  there  is  a 
weaker  thermocline  in  winter  than  in  summer.  In  the 
south  the  equatorial  climate  and  summer  inflow  from 
the  shallow  Java  Sea  cause  the  water  mass  to  be  fresher, 
warmer,  and  with  a  slightly  deeper  thermocline.  The 
model  SCS  summer  thermocline  depth  is  in  good  agree¬ 
ment  with  Ma  (1998),  at  30  and  50  m,  as  opposed  to 
the  model  winter  depth,  which  is  much  shallower.  The 
discrepancy  in  winter  may  be  caused  by  using  restoring 
type  surface  fluxes. 

c.  Volume  transport 

The  volume  transport  streamfunctions  for  summer 
and  winter  are  shown  in  Figs.  11a  and  1  lb.  For  the  SCS 
deep  basin,  the  barotropic  flow  reveals  an  anticyclonic 
(cyclonic)  gyre  in  the  summer  (winter).  For  the  SCS 
continental  shelf/slope  regions,  the  volume  transport  has 
an  evident  seasonal  variation.  Winter  volume  transport 
is  strongest  along  the  north  and  west  slope,  with  sig¬ 
nificant  eddy  and  meander  activity,  and  across  the  south¬ 
ern  slope  of  the  deep  basin.  Average  winter  VCJ  trans¬ 
port  is  about  10.6  Sv.  Summer  shows  a  weaker  VCJ 
transport  (5.5  Sv)  with  similar  transport  along  the  south¬ 
ern  slope  of  the  deep  basin,  again  with  noticeable  eddy 
activity  (Table  2). 

d.  Natuna  Island  eddy 

During  winter,  the  confluence  of  the  VCJ  and  the 
offshore  flow  from  Borneo  in  the  Sunda  shelf  generates 
a  mesoscale  cyclonic  eddy  (Fig.  6b),  here  termed  the 


Natuna  Island  eddy  (NIE),  consistent  with  observations 
of  late  winter  surface  cyclonic  flow  near  Natuna  Island 
described  by  Wyrtki  (Fig.  3b).  This  eddy  shows  little 
variability  during  the  winter,  with  maximum  swirl  ve¬ 
locity  of  0.6m  s_1  at  the  peak  of  the  winter  monsoon 
but  no  variation  in  position.  During  the  summer,  NIE 
becomes  a  weak  anticyclonic  eddy  (Fig.  6a).  As  shown 
in  Fig.  12,  the  winter  NIE  has  an  average  core  velocity 
of  0.45-0.5  m  s-1.  The  east-west  slanted  eddy  core  is 
deeper  on  the  west  side  (depth  —40  m)  of  the  eddy  than 
on  the  east  side  (depth  —20  m).  The  vertical  shear  of 
horizontal  velocity  promotes  baroclinic  instability. 

4.  Driving  mechanisms 

We  analyzed  the  results  of  the  four  experiments  to 
identify  the  driving  mechanisms  for  the  SCS  circulation. 

a.  Effects  of  nonlinearity  ( run  1  —  run  2 ) 

In  the  first  sensitivity  study  the  nonlinear  advection 
terms  were  removed  from  the  dynamic  equations.  Oth¬ 
erwise,  the  same  parameters  as  for  the  control  run  were 
used.  The  near-surface  vector  velocities  for  the  summer 
(Fig.  13a)  and  winter  (Fig.  13b)  illustrate  the  similarity 
to  the  control  run  results  (Fig.  6),  which  indicates  that 
the  nonlinearity  does  not  change  the  general  circulation 
pattern.  However,  our  computation  shows  that  the  non¬ 
linearity  causes  a  noticeable  change  in  the  volume  trans¬ 
port  and  the  strength  of  the  western  boundary  currents. 

Figures  14a  and  15a  are  plots  of  the  difference  in  the 
volume  transport  streamfunction  (A  if/)  for  the  control 
run  minus  the  linear  run  and  represents  nonlinear  effects 
for  summer  and  winter,  respectively.  The  biggest  A  ijj 
during  the  summer  (winter)  season  is  a  deep-basin  dou- 
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Fig.  10.  Thermohaline  feature  simulated  by  the  control  run  at  the  13°N  zonal  cross  section:  (a)  summer  temperature  (°C),  (b)  winter  temperature  (°C),  (c)  summer  salinity  (psu),  (d)  winter 
salinity  (psu),  and  at  113°E  latitudinal  cross  section:  (e)  summer  temperature  (°C),  (f)  winter  temperature  (°C),  (g)  summer  salinity  (psu),  and  (h)  winter  salinity  (psu).  Contour  interval  for 
temperature  is  2°C  and  for  salinity  0.2  psu. 
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Fig.  11.  Mean  total  volume  transport  streamfunction  (Sv)  for  control  run  during  (a)  summer  and  (b)  winter.  Contour  interval  is 

5  Sv.  The  shaded  regions  indicate  negative  values. 


ble-gyre  feature:  an  an ticy clonic  gyre  with  a  minimum 
value  of  -4  Sv  (-6  Sv)  west  of  Luzon  (16°N,  116°E) 
and  a  cyclonic  gyre  with  a  maximum  value  of  8  Sv  (12 
Sv)  west  of  Liyue  Bank  (1 1°N,  1 15°E).  The  double-gyre 
structure  in  winter  indicates  (i)  a  westward  cross-basin 
transport  (15  Sv)  near  13°N  (between  the  two  gyres), 
(ii)  increase  (decrease)  of  cyclonic  (anticyclonic)  trans¬ 
port  in  the  southwestern  basin  (south  of  13°N),  and  (iii) 
increase  (decrease)  of  anticyclonic  (cyclonic)  transport 
in  the  northeastern  basin  (north  of  13°N).  Furthermore, 
two  nearly  parallel  zero  curves  of  A  i/j  show  up  near  the 
Vietnam  coast  during  the  summer  (Fig.  14a).  One  curve 
is  close  to  the  coast,  and  the  other  curve  is  offshore. 
This  shows  that  A  i/j  =  0  on  the  entire  Vietnam  shelf, 
resulting  in  no  change  in  the  northward  VCJ  transport 
in  the  summer.  However,  during  the  winter  such  a  fea¬ 
ture  (Ai/j  =  0)  appears  only  north  of  13°N,  and  some 
southward  volume  transport  increment  (Ai/j  =  4  Sv) 
shows  up  near  the  Vietnam  coast  south  of  13°N.  After 


Table  2.  Volume  transport  (Sv)  of  individual  features  simulated 
by  the  control  run. 


Feature 

Transport 

Summer  Vietnam  coastal  jet 

5.5 

Winter  Vietnam  coastal  jet 

10.6 

Natuna  Island  eddy 

5.0 

averaging  A  i/j  over  the  Vietnam  coast  shelf  break,  we 
estimated  that  the  nonlinear  advection  brings  a  3.1  Sv 
of  the  VCJ. 

In  contrast  to  the  situation  at  13°N,  nonlinear  effects 
in  summer  are  important  at  6°N.  A  cross  section  of  the 
latitudinal  velocity  difference  during  the  summer  at  6°N 
shows  that  nonlinear  advection  augments  the  coastal 
core  current  at  the  western  boundary  by  about  0.25  m 
s_1  and  is  responsible  for  offshore  eddy  activity  (Fig. 
16a).  Such  a  nonlinear  effect  was  also  found  in  the 
surface  velocity  vector  difference  (AV)  field.  For  both 
seasons,  the  nonlinear  effect  on  circulation  is  very  weak 
in  the  deep  SCS  basin  and  quite  evident  near  the  west 
coast  and  southern  shelf  regions  (Figs.  17a  and  18a). 
The  nonlinear  effect  is  essential  for  the  NIE  formation 
near  110°E,  6°N  (Fig.  18a). 

b.  Wind-induced  circulation  ( run  1  —  run  3) 

The  difference  between  Run  1  and  Run  3  shows  the 
monsoon  wind  effects  on  the  SCS  surface  circulation 
and  SST  for  summer  and  winter  respectively.  The  fea¬ 
ture  of  the  summer  wind  effect  (Fig.  17b)  on  circulation 
and  volume  transport  is  in  contrast  to  the  winter  wind 
forcing  (Fig.  18b).  The  winter  wind  forcing  causes  the 
cyclonic  gyre  northeast  Natuna  Island  (the  NIE)  with 
an  average  0. 3-0.4  ms-1  rotational  velocity,  the  cross- 
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Fig.  12.  Mean  latitudinal  velocity  for  control  run  at  the  6°N  cross  section  during  (a)  summer  and 
(b)  winter.  Contour  interval  is  0.1  m  s-1. 
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Fig.  13.  Mean  surface  circulation  for  linear  run  during  (a)  summer  and  (b)  winter. 
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Fig.  14.  Summer  volume  transport  anomaly  (Sv)  caused  by  effects  of  (a)  nonlinearity,  (b)  winds,  and  (c)  boundary  forcing.  Contour 

interval  is  2  Sv. 


basin  flow  and  connecting  loop  across  the  Sunda  shelf 
around  Natuna  Island.  The  winter  VCJ  flows  with  0.3- 
0.4  m  s_1  velocity  (Fig.  18b),  corresponding  to  a  south¬ 
ward  volume  transport  of  4.6  Sv  (Fig.  15b).  The  China 
slope  flow  is  actually  slightly  more  than  in  the  control 
run.  When  winds  are  removed  from  the  forcing,  the 
model  results  show  considerable  change  in  the  structure 
of  the  summer  circulation  and  a  virtual  disappearance 
of  features  in  the  winter  circulation  (note  similarities  in 
Figs.  6b  and  18b);  surface  velocities  are  reduced  by 
0.20-0.3  ms”1. 


The  summer  winds  drive  a  mesoscale  anticyclonic 
eddy  near  a  shallow  depression  in  the  Sunda  shelf  and 
Natuna  Island,  with  rotation  speeds  of  0.15-0.2  m  s_1 
(Fig.  16b.)  Winds  also  induce  the  loop  current  joining 
the  Sunda  shelf  flow  and  the  VCJ,  bifurcation  of  the 
coastal  current,  and  cross-basin  flow  (Fig.  17b.)  The 
corresponding  wind  effect  on  summer  volume  transport 
is  largely  limited  to  the  northern  SCS  (Fig.  14b)  with 
a  6-8  Sv  maximum  transport  along  the  China  slope  and 
a  4-6  Sv  eastward  transport  centered  on  13°N,  a  region 
of  significant  nonlinear  effect. 


Fig.  15.  Winter  volume  transport  anomaly  (Sv)  caused  by  effects  of  (a)  nonlinearity,  (b)  winds,  and  (c)  boundary  forcing.  Contour  inter¬ 
val  is  2  Sv. 
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Fig.  16.  Summer  latitudinal  velocity  anomaly  (m  s_1)  at  the  6°N  cross  section  and  the  west  coast  zoomed  view,  caused  by  effects  of 
(a)  nonlinearity,  (b)  winds,  and  (c)  boundary  forcing.  Positive  (negative)  values  indicate  northward  (southward)  flow;  contour  interval 
is  0.1  m  s-1 . 


(a)  (b)  (c) 

Summer  (Cntrl-Linear)  Summer  (Cntrl-No  Wind)  Summer  (Cntrl-Cisd  Bdys) 


Longitude  (E)  Longitude  (E)  Longitude  (E) 

Fig.  17.  Summer  surface  current  anomaly  caused  by  effects  of  (a)  nonlinearity,  (b)  winds,  and  (c)  boundary  forcing. 
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Fig.  18.  As  in  Fig.  17  but  for  the  winter  surface  current  anomaly. 


Summer  wind  reduces  SST  slightly  (1°-2°C)  from 
up  welling  off  central  and  northern  Vietnam  and  the  east 
coast  of  Hainan  (Fig.  19a).  Winter  wind  decreases  sur¬ 
face  temperature  to  a  lesser  extent  on  the  east  side  of 
the  basin  along  the  Palawan  Trough  (Fig.  19b). 

Summer  winds  shift  the  average  surface  elevation  to 
the  north,  increasing  the  elevation  5-10  cm,  while  else¬ 
where  decreasing  the  height  by  the  same  magnitude 
(Fig.  20a).  The  winter  wind  increases  the  surface  ele¬ 
vation  over  the  continental  shelves  by  0.15-0.18  m,  and 
decreases  the  elevation  over  the  deep  basin  (Fig.  21a). 


Summer  SST  fCntrl-No  Wind) 


c.  Inflow/outflow  induced  circulation  (run  1  —  run  4) 

The  third  sensitivity  study  used  the  control  run  equa¬ 
tions  and  forcing  but  closed  all  open  lateral  boundaries, 
preventing  transport  of  mass,  heat,  or  salinity  through 
the  Luzon  Strait,  Taiwan  Strait,  or  Gasper  and  Karimata 
Straits.  With  no  inflow  or  outflow  the  summer  anticy- 
clonic  gyre  and  winter  cyclonic  gyre  are  more  pro¬ 
nounced.  Increased  recirculation  generally  leads  to 
greater  horizontal  and  vertical  variability  of  the  current 
structure. 


Winter  SST  (Cnfrl-No  Wind) 


Fig.  19.  Sea  surface  temperature  anomaly  caused  by  wind  effect.  Contour  interval  is  0.1°C;  the 
shaded  regions  indicate  negative  values. 


November  1999 


NOTES  AND  CORRESPONDENCE 


2987 


Longitude  (E)  Longitude  (E) 

Fig.  20.  Summer  surface  elevation  anomaly  (m)  caused  by  effects  of  (a)  winds,  and  (b) 
boundary  forcing.  Contour  interval  is  0.05  m. 


The  results  for  summer  lateral  boundary  forcing  retain 
more  of  the  anticyclonic  nature  of  the  circulation  than 
do  those  from  wind  forcing;  see  volume  transport 
streamfunction  difference  A i/s  (Figs.  14b,c).  The  flow 
generally  hugs  the  slope.  The  largest  gradient  of  Ai/jis 
located  off  the  south  Vietnam  coast  (Fig.  14c),  causing 
the  strong  western  boundary  currents.  The  southern  in¬ 
flow  splits  south  of  Natuna  Island  and  fans  across  the 
Sunda  shelf,  forming  a  small  anticylonic  eddy  collo¬ 


(a) 

Winter  Elv  (Cntrl-No  Wind) 


Longitude  (E) 


cated  with  that  forced  by  the  summer  wind.  The  VCJ 
transports  5.5  Sv  and  is  accompanied  by  a  countercur¬ 
rent  nearby.  At  6°N  (Fig.  16c),  however,  notice  that  the 
average  northward  velocity  throughout  most  of  the  shal¬ 
low  shelf  region  shows  no  discernible  pattern  of  change, 
while  the  shelf  current  along  the  100  m  isobath  de¬ 
creases  by  0.05-0.1  m  s_1.  The  summer  circulation  pat¬ 
terns  change  dramatically  when  the  lateral  boundary 
transport  is  removed.  This  can  be  seen  from  the  surface 


(b) 

Winter  Elv  (Cntrl-CIsd  Bdys) 


Longitude  (E) 


Fig.  21.  Winter  surface  elevation  anomaly  (m)  caused  by  effects  of  (a)  winds,  and  (b)  boundary 
forcing.  Contour  interval  is  0.05  m. 
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velocity  vector  difference  AV  (Fig.  17c).  There  is  still 
western  intensification  of  the  current  along  the  coast  of 
Malaysia,  and  it  still  joins  the  flow  out  of  the  Gulf  of 
Thailand  to  contribute  to  an  intensification  of  current 
off  southern  Vietnam. 

Winter  closed  boundary  circulation  patterns,  on  the 
other  hand,  show  less  difference  in  structure  from  the 
control  run  but  more  variability  in  magnitude.  There  is 
an  average  decrease  of  0. 1-0.4  m  s_1  in  current  speed 
of  the  winter  VCJ,  which  equates  to  4-6  Sv  southward 
volume  transport  (Fig.  15c).  The  Kuroshio  intrusion  and 
inflow  through  Taiwan  Strait  are  obviously  supple¬ 
mented  by  recirculation  flow  from  along  the  coast  of 
Luzon  Island.  Near  Natuna  Island  it  is  noteworthy  that 
the  structure  of  cross-basin  circulation  and  current  flow 
away  from  the  Borneo  coast  are  unchanged.  Similarly, 
it  is  apparent  that  the  spatial  extent  and  shape  of  the 
gyre  northeast  of  Natuna  Island  is  unchanged  by  closing 
the  boundary  flow.  In  cross  section,  however,  it  can  be 
seen  that  the  NIE  does  lose  some  of  its  velocity — there 
is  a  0. 1-0.3  m  s-1  decrease  in  the  average  core  velocity 
on  the  western  side  at  6°N  in  this  run  corresponding  to 
a  volume  transport  of  3  Sv,  associated  with  the  decrease 
in  velocity  of  the  VCJ.  The  difference  in  velocity  is 
much  smaller  on  the  eastern  side,  suggesting  that  some 
other  effect  dominates  the  current  structure  in  that  re¬ 
gion. 

The  summer  lateral  boundary  transport  decreases  the 
average  surface  elevation  over  the  China-Vietnam  con¬ 
tinental  shelf  by  0.05-0.15  m,  while  elsewhere  increases 
the  height  with  a  maximum  accretion  of  0. 2-0.4  m  near 
the  Karimata  Strait  (Fig.  20b).  The  winter  lateral  bound¬ 
ary  transport  increases  the  average  surface  elevation 
over  the  China-Vietnam  continental  shelf  by  0.05-0.2 
m,  while  elsewhere  decreases  the  height  with  a  maxi¬ 
mum  reduction  of  0.25  m  near  the  Karimata  Strait  (Fig. 
21b). 

5.  Conclusions 

1)  The  SCS  circulation  and  thermohaline  structure 
was  simulated  in  this  study  by  the  POM  model  under 
the  climatological  forcing.  During  the  summer  (winter) 
monsoon  period  the  SCS  surface  circulation  is  generally 
anticyclonic  (cyclonic)  with  a  strong  western  boundary 
current  (width  around  100  km) — the  Vietnam  coastal 
jet.  This  jet  has  a  strong  seasonal  variation  and  flows 
northward  (southward)  during  summer  (winter)  with  a 
mean  maximum  speed  of  0.5  m  s-1  (0.95  m  s_1),  with 
a  mean  volume  transport  of5.5Sv(10.6Sv)  extending 
to  a  depth  of  around  200  m  (500  m).  During  summer, 
the  western  boundary  current  splits  and  partially  leaves 
the  coast;  the  bifurcation  point  is  at  14°N  in  May  and 
shifts  south  to  10°N  in  July.  The  cross-basin  zonal  cur¬ 
rent  reaches  12°N  and  a  core  speed  of  0.6  m  s_1  by  the 
end  of  the  summer.  The  coastal  branch  continues  north, 
then  east  at  Hainan  Island.  Average  summer  sea  surface 
elevation  varies  from  —0.1  to  0.1  m  with  maximum 


values  at  northeastern  part  near  Luzon.  Sea  surface 
heights  for  the  winter  show  a  0.2-m  depression  over  the 
deep  basin  with  maximum  values  in  the  southwestern 
part  over  the  Sunda  shelf  and  the  Gulf  of  Thailand.  This 
suggests  that  the  southwest  monsoon  wind  in  boreal 
summer  piles  up  the  water  at  the  northeastern  part  near 
Luzon  and  the  northeast  monsoon  wind  in  boreal  winter 
piles  up  the  water  at  the  southwestern  part  over  the 
Sunda  shelf  and  the  Gulf  of  Thailand.  The  POM  model 
also  successfully  simulates  a  mesoscale  eddy  in  the  Sun¬ 
da  shelf,  namely  the  Natuna  Island  eddy.  This  eddy  is 
cyclonic  (anticyclonic)  with  maximum  swirl  velocity  of 
0.6  m  s-1  at  the  peak  of  the  winter  (summer)  monsoon. 

Isotherms  and  isohalines  for  summer  and  winter  are 
nearly  horizontal  from  east  to  west  except  at  the  coastal 
regions.  Coastal  up  welling  and  downwelling  are  also 
simulated:  localized  lifting  (descending)  of  the  iso¬ 
therms  and  isohalines  during  summer  (winter)  at  the 
west  boundary.  Both  coastal  up  welling  and  downwell¬ 
ing  (causing  horizontal  thermohaline  gradient),  com¬ 
bined  with  the  high  velocity  shear  across  the  coastal  jet, 
results  in  baroclinic  instability.  This  mechanism  may 
contribute  to  the  summer  jet  bifurcation.  In  general  mod¬ 
el  thermohaline  structure  is  consistent  with  the  two  SCS 
water  masses  described  by  Wyrtki  (1961.)  Over  the 
southern  basin  there  is  a  general  lifting  of  isotherms  of 
40-50  m  from  winter  to  summer  above  200  m.  In  the 
northern  SCS,  near-surface  waters  (20-50  m)  are  influ¬ 
enced  by  the  winter  inflow  of  North  Pacific  Kuroshio 
water. 

2)  The  model  wind  effects  on  the  SCS  are  more  per¬ 
vasive  than  those  of  nonlinear  dynamic  effects.  Similar 
to  nonlinear  effects  there  are  areas  of  anticy Ionic  volume 
transport  midbasin  in  the  summer  and  winter,  although 
the  midbasin  surface  currents  reverse.  Summer  wind 
effects  on  transport  are  generally  divided  along  the  axis 
of  the  monsoon  wind,  SW-NE,  with  anticyclonicity  and 
positive  elevation  in  the  southeast  and  cyclonic  transport 
and  negative  elevation  in  the  northwest.  The  winter  wind 
effect  is  largely  cyclonic,  including  a  center  (with  a 
strong  surface  current  signature)  for  the  NIE,  with  a 
matching  steeper  depression  in  surface  elevation. 

3)  The  model  boundary  forcing  effect  on  volume 
transport  is  most  easily  distinguished  as  anticyclonic  in 
the  summer  and  cyclonic  in  the  winter.  The  stronger 
transport  centers  of  both  seasons  have  strong  surface 
current  signatures  including  the  NIE.  Topographically 
linked  eddy  flow  and  transport  is  evident  in  the  summer, 
while  winter  flow  narrowly  follows  the  200-m  isobath 
until  it  encounters  the  Sunda  shelf  and  Natuna  Island. 

In  the  summer,  SCS  volume  transport  is  mostly  due 
to  boundary  forcing.  Surface  movement  is  marginally 
dominated  by  wind  forcing,  when  the  boundary-forcing 
surface  effect  (away  from  the  southern  boundary)  co¬ 
alesces  with  associated  volume  transport.  Winds  drive 
the  offshore  summer  bifurcation  of  the  VCJ.  The  anti- 
cylonic  gyre  south  of  Hainan  is  a  surface  phenomena 


November  1999 


NOTES  AND  CORRESPONDENCE 


2989 


driven  by  boundary  flow  interaction  with  topography 
and  nonlinear  dynamics. 

In  the  winter,  surface  movement  due  to  boundary 
forcing  is  largely  limited  to  the  perimeter  of  the  SCS 
basin,  specifically  a  narrow  swath  over  the  200-m  iso¬ 
bath,  away  from  the  southern  boundary.  Surface  flow 
and  volume  transport  are  dominated  by  the  wind  effect 
except  at  the  NIE,  a  mesoscale  feature  that  must  be 
qualified  by  nonlinear  dynamics.  Also  to  be  qualified 
by  nonlinear  dynamics  is  cross-basin  transport  for  both 
seasons. 

4)  Future  studies  should  concentrate  on  less  simplis¬ 
tic  scenarios.  Realistic  lateral  transport  (especially  at  the 
Mindoro  and  Karimata  Straits)  and  surface  heat  and  salt 
fluxes  should  be  included,  and  the  use  of  extrapolated 
climatological  winds  needs  to  be  upgraded  to  incor¬ 
porate  synoptic  winds  to  improve  realism.  Finally,  the 
assumption  of  quasi  linearity  that  allowed  us  to  use 
simple  difference  to  quantify  the  effect  of  external  forc¬ 
ing  needs  to  be  rigorously  tested.  It  is  very  important 
to  develop  a  thorough  methodology  to  perform  sensi¬ 
tivity  studies  under  the  highly  nonlinear  conditions  that 
may  exist  in  the  littoral  environment. 
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